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PREFACE 


The  Structures  and  Materials  Panel  of  AGARD  has  sponsored  a  number  of  activities 
concerned  with  the  behaviour  of  materials  at  elevated  temperatures  in  gas  turbines.  These 
have  covered  a  number  of  aspects  from  laboratory  investigations  to  the  practical  difficulties 
of  service  maintenance  of  gas  turbine  engine  parts.  This  is  demonstrated  in  the  Specialists’ 
Meetings  held  by  the  SMP,  namely: - 

“Characterisation  of  Low  Cycle  High  Temperature  Fatigue  by  the  Strain  Range 

Partitioning  Method”  (Aalborg  1978,  AGARD  CP-243) 

“Ceramics  for  Gas  Turbine  Applications”  (Porz-Wahn  1979,  AGARD  CP-276) 

“Maintenance  in  Service  of  High  Temperature  Parts”  (Noordwijkerhout  1981 , 

AGARD  CP-317). 

In  the  final  discussion  of  the  last  Specialists’  Meeting  there  was  a  clear  indication  that 
more  materials  performance  data  under  engine  conditions  was  needed  by  the  user.  There  was 
also  a  requirement  for  some  basis  of  damage  prediction  upon  which  to  build  an  approach  of 
damage  tolerance  design  and  use  of  engine  components.  In  response  to  this  the  SMP  approved 
a  follow-on  activity  entitled  “Damage  Tolerance  Concepts  for  Critical  Engine  Components”. 

One  of  the  early  activities  of  the  group  involved  was  to  receive  a  presentation  from 
Dr-Ir.  H.P.  van  Leeuwen  of  NLR  on  “The  Application  of  Fracture  Mechanics  to  the  Growth 
of  Creep  Cracks”  in  London,  Spring  1983.  I  am  pleased  to  see  this  significant  review  published 
as  an  AGARD  Report. 


D.A.FANNER 

Chairman  -  Sub-Committee  on 
Damage  Tolerance  Concepts  for  Critical 
Engine  Components 
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by 

Dr.Ir.  H.P.  van  Leeuwen 
National  Aerospace  Laboratory  NLR 
P.O.  Box  90502,  1006  BM  AMSTERDAM 
The  Netherlands 


SUMMARY 

— The  literature  concerning  the  application  of  fracture  mechanics  to  the  growth  of  creep  cracks  is 
discussed. 

As  an  introduction  it  is  shown  that  linear  elastic  fracture  mechanics  have  been  applied  successfully 
to  the  residual  strength  of  brittle  materials,  crack  growth  due  to  stress  corrosion  and  crack  growth  due  to 
fatigue. 

Next  the  methods  are  discussed  to  account  for  plasticity. 

Then  a  survey  is  made  of  the  attempts  to  correlate  the  crack  growth  rate  under  creep  conditions  with’ 
the  stress  intensity  factor,  the  net  section  stress,  the  reference  stress,  the  crack  opening  displacement 
rate,  the  contour  integral  J  and  the  contour  integral  J*. 

A  few  sections  are  devoted  to  alternative  methods  to  calculate  J  and  J* . 

A  survey  is  made  also  of  theoretical  predictions  of  the  crack  growth  rate  as  a  function  of  parameters 
that  are  considered  to  be  determining. 

The  report  ends  with  a  discussion  and  a  few  conclusions,  'rr" 


LIST  OF  ABBREVIATIONS 

ASTM  American  Society  for  Testing  and  Materials 

CCP  Centre  Cracked  Plate 

CCT  Centre  Cracked  Tension 

CDCB  Contoured  Double  Cantilever  Beam 

CN  Centre  Notch 

CNC  Centre  Notch  Cylindrical 

CNP  Centre  Notched  Plate 

CNRS  Circumferentially  Notched  Round  Specimen 

CNT  Centre  Notch  Tension 

COD  Crack  Opening  Displacement 

CT  Compact  Tension 

CTOA  Crack  Tip  Opening  Angle 

CTOD  Crack  Tip  Opening  Displacement 

DEN  Double  Edge  Notch 

DENT  Double  Edge  Notch  Tension 

FEM  Finite  Element  Method 

LEFM  Linear  Elastic  Fracture  Mechanics 
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RNB  Round  Notched  Bar 
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t .  INTRODUCTION 

The  last  two  decades  have  seen  the  development  of  an  aircraft  design  philosophy  based  on  the 
principle  of  damage  tolerance.  This  means  that  the  fact  is  accepted  that  the  occurrence  of  cracks  or  other 
forms  of  damage  in  the  aircraft  structure  cannot  always  be  prevented. 

It  is  then  required  that  the  cracked  component  can  sustain  a  certain  fraction  of  its  design  load  without 
collapse  due  to  unstable  crack  growth.  It  is  required  also  that  in  service  the  crack  will  grow  so  slowly 
that  it  can  be  detected  by  periodic  inspections  before  it  reaches  a  catastrophic  size.  If  a  crack  is 
expected  in  a  location  which  is  inaccessible  for  inspection  it  will  have  to  grow  so  slowly  that  it  will  not 
become  catastrophic  during  the  entire  design  life  of  the  structure.  An  alternative  design  measure  is  to 
provide  a  neighbouring  component  that  will  take  the  load  from  a  failing  component  without  failing  itself. 
This  design  philosophy  requires  calculations  leading  to  prediction  of  crack  growth  rates  and  residual 
strength.  Fracture  mechanics  provide  the  tools  fo.r  such  calculations. 

Fracture  mechanics  have  been  applied  successfully  to  airframe  structures  where  crack  growth  is  mainly 
by  fatigue  at  moderate  or  normal  temperatures. 

In  recent  years  the  desire  has  grown  to  apply  principles  similar  to  damage  tolerance  to  aircraft  jet 
engines.  This  has  raised  the  question  whether  fracture  mechanics  can  be  applied  to  crack  growth  by  creep 
and  fatigue  in  engine  components  operating  at  high  temperatures . 

A  similar  problem  has  arisen  with  regard  to  stationary  high  temperature  installations  such  as 
furnaces,  chemical  reactors  and  boilers,  especially  those  in  electric  power  plants.  Such  installations  have 
been  designed  typically  for  a  finite  life  and  on  the  basis  of  creep  rupture  strength.  Many  of  such 

installations  have  now  reached  their  design  life  time.  For  economic  reasons  their  owners  wish  to  keep  these 

installations  in  service.  The  question  is  how  long  this  can  be  done  and  how  safety  can  be  guaranteed.  An 

important  aspect  is  that  in  some  of  these  installations  cracks  have  been  found  and  that  present  require¬ 

ments  do  not  allow  cracks. 

This  has  spurred  many  investigations  regarding  the  applicability  of  fracture  mechanics  to  crack  growth 
under  creep  conditions. 

Here  the  problem  is  that  conventional  fracture  mechanics  are  based  on  linear-elastic  stress-strain 
relations.  Because  the  used  structural  materials  show  a  certain  measure  of  plasticity  before  failure  one 

was  forced  to  account  for  plasticity  effects. 

Many  equations  in  linear  elastic  fracture  mechanics  (LEFM),  can  if  necessary  by  some  adaptation,  be  applied 
provided,  however,  the  region  where  plasticity  occurs  remains  small  relative  to  the  length  of  the  crack 
and  the  dimensions  of  the  specimen  or  the  component.  Plastic  fracture  mechanics,  for  cases  where  this 
condition  is  not  fulfilled,  are  still  under  development. 

In  the  creep  case  one  is  dealing  with  an  important  amount  of  time-dependent  plastic  deformation  which 
renders  the  applicability  of  conventional  fracture  mechanics  highly  doubtful. 

In  the  past  15  years  or  so  several  empirical  and  theoretical  studies  have  been  made  into  the  possibility 
of  applying  fracture  mechanics  to  the  growth  of  creep  cracks. 

Reviews  of  such  studies  have  been  published  earlier  by  Haigh  (1),  Van  Leeuwen  (2,  3),  Ohji  (4),  Ellison 
and  Harper  (5),  Ashby  and  Tomkins  (6),  Pilkington  (7),  Fu  ?8),  Radhakrishnan  and  McEvily  (9)  and  by 
Sadananda  and  Shahinian  (10). 

The  present  review  tries  to  be  more  complete  and  will  treat  some  more  recent  literature, 
i  As  an  introduction  the  application  will  be  treated  of  LEFM  to  the  residual  strength  of  brittle  materials, 

crack  growth  by  stress  corrosion  and  crack  growth  by  fatigue.  Next  a  review  will  be  given  of  methods  to 
;  account  for  plasticity,  largely  taken  from  an  article  by  Eftis  and  Liebowitz  (11).  Then  a  review  is 

,  presented  of  attempts  to  correlate  creep  crack  growth  rate  with  the  stress  intensity  factor  K,  the  net 

,  section  stress  o  ,  the  reference  stress  a  the  crack  (tip)  opening  displacement  rate  v  (or  5),  the 

•  contour  integral0  J  and  the  contour  integraf  J*.  The  report  ends  with  conclusions  and  recoranendat ions . 

I 

i  2.  LINEAR  ELASTIC  FRACTURE  MECHANICS 

.  LEFM  has  found  a  successful  application  in  the  description  and  the  prediction  of  several  types  of 

!  crack  growth.  A  condition  is  that  the  extent  of  plastic  deformation  will  remain  small.  This  implies  that 

the  plastic  zone  at  the  crack  tip  will  remain  small  and  will  remain  embedded  in  the  field  where  deformation 
is  elastic  only.  In  this  way  the  plastic  deformation  will  have  a  negligible  effect  on  the  stresses  in  the 
elastic  field  and  on  the  quantities  calculated  with  them. 

2.1  Residual  Strength 

■  If  LEFM  is  applied  to  a  residual  strength  test  on  a  cracked  specimen  it  is  assumed  that  unstable  crack 

'  growth  will  occur  if  the  stress  intensity  factor  reaches  a  critical  value  K  .  It  is  then  found  that  K  is 

not  really  a  material  constant  independent  of  specimen  shape  and  size  but  will  depend  heavily  on  specimen 

‘  thickness.  See  for  instance  figure  1  taken  from  Sullivan  and  Freed  (12). 

2.2  Stress  Corrosion  Cracking 

|  In  the  application  of  LEFM  to  SCC  it  is  found  that  the  crack  growth  rate  depends  on  K  in  a  very 

|  typical  way  (Fig.  2,  Speidel ,  13)  and  shows  a  region  I  where  depends  strongly  on  K,  a  region  II  where 

^  is  constant  (velocity  plateau)  and  a  region  III  where  ^  again  rises  rapidly  with  K.  In  practice 

the  relationship  between  da/d t  and  K  appears  to  be  influenced  extensively  by  such  conditions  as  material 
j  susceptibility,  relative  humidity,  halide  concentration,  pH-vslue,  electrode  potential,  temperature, 

,  viscosity  of  the  corrodent,  nature  of  the  corrodent  (e.g.  liquid  metal  v.s.  water). 

I 

I  2.3  Fatigue 

I 

The  application  of  LEFM  to  SCC  is  possible  because  there  stresses  usually  remain  small.  This  is 
'  usually  also  the  case  with  fatigue.  In  this  case  the  crack  growth  rate  da/dN  is  determined  by  the  variation 
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AK  of  K.  The  relationship  is  influenced  rather  heavily,  however,  by  frequency  and  environment  (Fig.  3, 
Speidel,  14). 

2.4  Conclusion 

Although  certain  relationships  have  been  established  between  fracture  mechanics  parameters  and 
various  types  of  cracking,  these  relationships  are  influenced  strongly  by  a  number  of  additonal  conditions 
or  parameters.  If,  therefore,  a  correlation  can  be  found  in  the  case  of  creep  cracking  one  must  be  prepared 
to  accept  that  also  in  this  case  a  number  of  additional  parameters  can  influence  the  correlation  to  a  large 
extent. 


3.  PLASTICITY  EFFECTS 

In  many  technical  materials  appreciable  plastic  deformation  can  occur  before  the  crack  becomes  un¬ 
stable.  If  the  plastic  zone  at  the  crack  tip  remains  small  and  remains  embedded  in  the  elastic  field 
certain  corrections  can  be  applied  to  the  equations  used  to  calculate  the  stress  intensity  factor  K. 

When  plastic  deformation  becomes  more  extensive,  the  K  factor  may  lose  its  validity  and  one  tries  to 
apply  other  parameters  that  can  be  expected  to  give  a  better  description  of  crack  tip  conditions.  Such 
parameters  are  for  instance 

-  the  crack  (tip)  opening  displacement,  C(T)OD 

-  the  J  integral 

-  the  net  se.ction  stress 

These  parameters  will  now  be  reviewed  because  they  provide  suggestions  for  the  treatment  of  creep  cracking 
by  fracture  mechanics  methods. 


3.1  Effective  Crack  Length 


The  value  of  the  K-factor  depends  essentially  on  the  square  root  of  the  crack  length  a.  In  the  case 

y  i 

lat 

„2 


of  plasticity  it  can  be  assumed  that  the  plastic  zone  at  the  crack  tip  is  a  circle  with  radius  r  where  r 
can  be  calculated  as  ' 


B§- 

°Y 


where  6  -  1/2ir  for  the  plane  stress  case  and  B  ”  1  /6r  for  the  plane  strain  case. 
Plasticity  is  then  acccounted  for  by  introduction  of  an  effective  crack  length 


a  +  r„ 


(3.1) 


(3.2) 


In  this  way  the  effect  plastic  deformation  has  on  the  elastic  stress  field  1b  accounted  for  in  a 
approximative  manner. 

3.2  Energy  Dissipation 

The  K-value  is  a  scaling  factor  for  the  elastic  stress  field  at  the  crack  tip  that  can  be  represented 


°jk 
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(3.3) 


1 2irr 


In  reality  this  is  the  first  term  of  a  series  expansion  with  increasing  powers  of  the  quantity 
"  -  1 

(r/a)2  with  n  ■  1,  2,  3,  ... 

Higher  values  of  (j  -  1)  can  be  neglected  if  r  remains  sufficiently  small  relative  to  a,  e.g.  for  r  <  0.05a. 
The  original  fracture  mechanics  are  based  on  energy  considerations  however. 

A  characteristic  quantity  is  the  strain  energy  release  rate 


G 

da 


(3.4) 


where  U  is  the  potential  energy. 

For  the  elastic  case  G  can  be  calculated  as 


■4 
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(3.5) 


Here  B  is  the  specimen  thickness,  C  the  compliance  v/P  wherein  v  is  the  load  line  displacement  and  P  the 
load. 

It  can  be  shown  that  for  the  elastic  case  G  is  proportional  to  K  . 

If  the  extent  of  plasticity  is  so  great  that  the  K-factor  loses  its  significance,  G  may  still  be  a  useful 
parameter  to  describe  the  material's  resistance  to  cracking. 

In  thin  sheet,  where  the  plastic  zone  may  be  larger  than  the  sheet  thickness,  stable  crack  growth  will 
occur  when  the  gross  section  stress  a  is  raised  monotonically.  The  resistance  to  stable  crack  growth  can  be 
depicted  by  the  R-curve  (Fig.  4). 

If  on  a  specimen  with  a  precrack  of  length  a^  the  stress  o  is  raised,  G  will  rise  in  proportion  with  a  and 
a.  Every  stress  increase  Ao  will  result  in  a  crack  length  increase  Aa.  The  locus  of  the  points  on  the  G-lines 
corresponding  to  a  certain  stress  o.  and  the  resulting  crack  length  a.  forms  the  R-curve. 

Fracture  instability  will  occur  whe/l  the  C-line  is  tangent  to  the  R-curve  because  this  implies  that  at 
constant  stress  a  more  energy  is  released  than  is  needed  for  crack  growth  since 
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Related  to  this  approach  is  the  one  by  Eftis  and  Liebowitz  (11)  who  calculate  the  energy  release  when  the 
load-displacement  or  P-v  diagram  is  non-linear  due  to  crack  growth  and  plasticity  (Fig.  5). 

They  represent  the  non-linear  curve  by  the  equation 

v  "  5  +  k(5)D  (3‘7) 

where  M  is  the  initial  rigidity. 

Through  application  of 

G--$E  (3.8) 

da 

they  arrive  at  the  effective  energy  release  rate  G  as 
n,  t  /^\n-1i 


Unstable  crack  growth  will  occur  for  a  critical  energy  release  rate  Gc> 

3.3  J-integral 

Related  to  the  approach  based  on  G  is  the  one  based  on  the  contour  integral  J,  proposed  by  (Jim)  Rice 
(Fig.  6) 


W  is  the  strain  energy  density 

e 

mn 

wV  °jkd£jk  (3'n) 

T.  is  a.  n.  ,  the  traction  normal  to  the  contour  T  and  u.  is  the  displacement  vector. 

J-^can  be  determined  for  the  elastic/plastic  case  through^a  numerical  finite  element  calculation  or 
analytically  e.g.  using  a  polar  coordinate  system  around  the  crack  tip. 

It  has  been  show  that  H  is  path  iundependent  provided  the  contour  stays  away  from  the  crack  tip  sufficiently. 
Proof  has  been  given  however,  using  the  so  called  deformation  theory  of  plasticity  which  is  valid  only  for 
monotonically  rising  stresses. 

It  has  been  shown  that  J  is  also  a  measure  of  the  energy  release  rate 


For  the  elastic  case  one  finds  J  »  G. 

J  can  be  determined  empirically  from  a  non-linear  diagram  of  load  P  versus  load  line  displacement  v, 
see  figure  7,  (Landes  and  Begley,  15).  It  is  postulated  that  further  crack  growth  will  occur  for  a 

critical  J  value  J  . 

c 

3.4  Crack  (Tip)  Opening  Displacement 

A  seemingly  different  approach  is  the  one  based  on  crack  (tip)  opening  displacement. 

The  COD  can  be  measured  in  the  load  line,  or  at  other  distances  from  the  crack  tip,  using  e.g.  clip  gauges. 
By  calculation  and/or  calibration  the  CTOD  can  be  derived. 

It  can  be  calculated  using  the  plastic  crack  tip  model  of  Dugdale 


8o„a  /  \ 
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From  a  quasi-elastic  calculation,  using  a  circular  plastic  zone  around  the  tip  of  the  effective  crack 
(a  +  r^),  6  is  found  as  the  displacement  at  the  real  crack  tip. 

The  result  is 


where  a  ■  4/ir/3  for  plane  strain  and  a  •  4/tt  for  plane  stress. 
This  shows  the  relation  to  the  energy  dissipation  approach  since 


a  («T  '  6) 


3.5  Stable  Versus  Unstable  Crack  Growth 


When  plasticity  plays  a  role  appreciable  stable  crack  growth  can  occur  prior  to  failure. 

To  investigate  what  parameters  would  govern  or  describe  this  process  finite  element  calculations  have  been 
made  in  which  the  behaviour  of  v,  the  load  line  displacement,  was  forced  to  vary  with  a  in  the  manner  as 
had  been  determined  experimentally. 

De  Koning  (16)  found  that  re-initiation  of  the  crack  occurred  at  critical  values  Jc  or  6  of  J  and  6  but 
that  during  stable  crack  growth  the  crack  tip  opening  angle  (CTOA)  remained  virtually  constant. 

Shih  et  al.  (17)  found  that  dJ/da  remained  constant  over  a  short  stretch.  The  so  called  nominal  crack  tip 
opening  angle  d£/da  and  the  real  CTOA  remained  constant  over  a  large  span  of  a-values. 

They  concluded  that  crack  re-initiation  could  be  predicted  on  the  basis  of  J  or  6  but  that  for 
continued  crack  growth  CTOA  or  d6/da  had  to  have  a  certain  value. 

Gudas  et  al.  (18)  show  results  indicating  that  dJ/da  remains  fairly  constant  for  0.5  —  0.8. 

2  w 

They  conclude  that  Paris's  tearing  modulus  (dJ/da) . (E/oq  )  can  be  a  good  measure  for  the  onset  of 

instability. 
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3.6  Net  Section  Stress 

When  yielding  occurs  over  the  larger  part  of  the  reduced  section  it  is  expected  that  crack  instability 

will  occur  for  a  critical  value  of  the  net  section  stress  a  . . 

net 

4.  THE  GROWTH  OF  CREEP  CRACKS 

If  a  metal  specimen,  containing  a  crack  due  to  prior  fatigue  stressing  for  instance,  is  loaded  e.g. 
in  tension  or  in  bending  at  elevated  temperature,  the  crack  will  continue  to  grow,  perhaps  after  some 
incubation  period.  By  measuring  and  recording  the  crack  length,  using  the  change  in  electric  resistance 
or  in  compliance,  a  crack  growth  curve  as  depicted  in  figure  8  can  be  determined.  In  general  the  crack 
growth  curve  will  strongly  resemble  the  conventional  creep  elongation  curve  with  its  typical  three  regions 
with  fast,  slow  and  fast  growth  rates. 

In  the  case  of  the  creep  crack  growth  curve  the  primary  region  may  be  small  or  absent,  the  slope  in  the 
secondary  region  may  be  shallow  or  steep,  but  the  basic  form  will  always  be  there. 

We  are  dealing  again  with  stable  crack  growth  followed  by  unstable  fracture.  The  phenomenon  is  determined 
more  strongly  by  plasticity,  however,  as  can  be  seen  from  the  following  sunmary: 

-  Brittle  materials,  heavy  plate:  no  stable  crack  growth 

-  Ductile  materials,  thin  sheet:  stable  crack  growth  at  monotonically  increasing  load. 

-  Creeping  materials:  stable  crack  growth  at  constant  load. 

It  need  not  be  expected  that  final  failure  in  the  third  case  will  require  a  different  analytical  treatment 
than  that  in  the  second  case.  It  can  be  expected  however  that  stable  crack  growth  in  the  third  case  will 
have  to  be  treated  differently  than  that  in  the  second  case,  as  much  as  stable  crack  growth  in  the  second 
case  had  to  be  treated  differently  from  unstable  crack  growth  in  the  first  case. 

There  is  a  certain  resemblance  with 

-  fatigue:  stable  crack  growth  at  constant  cyclic  stress 

-  stress  corrosion:  stable  crack  growth  at  constant  static  stress. 

5.  CORRELATION  OF  THE  CREEP  CRACK  GROWTH  RATE 
WITH  THE  STRESS  INTENSITY  FACTOR,  THE  NET 
SECTION  STRESS  AND  THE  REFERENCE  STRESS 


As  was  explained  before,  the  applicability  of  the  K-factor  is  dependent  on  some  conditions,  the  most 
important  perhaps  that  the  extent  of  plasticity  should  be  relatively  small. 

It  is  somewhat  naive  may  be  to  suppose  that  K  will  be  applicable  to  a  crack  in  a  creep  specimen  where  time 
dependent  plastic  deformation  occurs  over  the  greater  part  of  the  net  section  and  perhaps  also  of  the 
gross  section  away  from  the  crack.  It  will  not  be  very  difficult  to  calculate  a  X-value  in  some  way  and 
to  make  a  diagram  of  da/dt  versus  K.  The  question  is  whether  this  would  be  meaningful.  For  that  it  is 
necessary  that  under  otherwise  identical  conditions,  for  different  specimen  shapes  and  sizes,  for  different 
crack  lengths  and  gross  section  stresses,  the  same  da/dt  will  be  found  if  K  is  the  same.  In  practice  this 
results  in  many  disappointments.  t 

If  a  diagram  of  secondary  creep  rate  versus  K  is  made  nevertheless,  a  curve  as  shown  in  figure  9  may  result. 
This  is  again  a  sigmoidal  curve.  One  distinguishes  a  region  I  where  da/dt  depends  strongly  on  K,  a  region  II 
where  da/dt  depends  weakly  on  K,  and  a  region  III  where  da/dt  again  rises  sharply  with  K. 

In  some  cases  the  primary  branch  may  be  absent  and  the  secondary  branch  may  be  horizontal.  In  other  cases 
the  secondary  branch  may  be  nearly  as  steep  as  the  primary  and  the  tertiary.  Some  investigators  show  only 
one  branch  and  than  it  is  difficult  to  judge  which  of  the  three  it  is. 

Note  that  the  crack  velocity  diagram  for  creep  has  some  resemblance  to  those  for  stress  corrosion  and 
fatigue. 

If  plasticity  is  extensive  the  correlation  of  da/dt  with  K  can  be  expected  to  be  bad.  One  can  then  attempt 
to  correlate  da/dt  with  net  section  stress  0ne(. .  For  a  cracked  specimen  loaded  partly  in  bending,  net 
section  stress  is  not  such  a  meaningful  thing  and  one  prefers  to  use  the  reference  stress  aref 
It  can  be  defined  as 


(5.0 


where  P„  is  the  load  at  the  onset  of  yielding  in  the  uncracked  specimen  and  Pc  is  the  load  at  the  onset  of 
yielding  in  the  cracked  specimen. 

Many  of  the  older  investigations  concentrated  on  the  question  whether  da/dt  would  correlate  well  with  K  and 
if  so  whether  it  would  correlate  better  with  K  than  with  o  or  o  .. 

A  typical  example  is  the  investigation  by  Ellison  and  WaltoS  (19;,  see  figure  10.  Single  edge  notch  flat 
specimens  have  been  loaded  axially  and  in  bending  at  a  series  of  initial  stresses.  It  is  seen  that  for  the 
axially  loaded  case  creep  crack  growth  rates  correlate  fairly  well  with  the  stress  intensity  factor. 

For  the  bending  case  there  is  a  large  spread  however,  and  plotting  da/dt  versus  K  does  not  reduce  the 
data  to  a  single  curve. 

Another  example  is  the  investigation  by  Neate  and  Siverns  (20),  figures  11  and  12.  A  0.5  Cr-0.5  Mo-0.25  V 
steel  has  been  given  two  different  heat  treatments,  one  involving  a  quench  making  the  steel  relatively 
btittle  in  order  t'  simulate  a  s  Id  heat  affected  some,  the  other  comprising  normalizing  and  tempering. 
Correlation  of  da.  .  with  X  i  asonably  good.  Somewhat  unexpectedly  the  correlation  with  o  is  rather 
lad,  especially  fo  ..  d-.sc:  t  material. 

A  third  example  invo’  AISi  3 1o  type  stainless  steel  (Fig.  13)  tested  by  Nicholson  and  Formby  (21). 

Here  the  correlation  with  7  is  better  than  with  K  as  one  would  expect  for  a  ductile  stainless  steel. 

As  was  shown  the  results  tofiltimes  do  live  up  to  expectations  and  sometimes  they  don’t.  Literature  in  which 
correlations  of  da/dt  with  K,  o  and  o  ,  were  compared  have  been  compiled  and  classified  as  to  the  result 

in  table  1.  It  is  seen  that  tl:e°results  are  rather  erratic  and  unpredictable.  Therefore  correlations  with 
either  o  or  K  should  not  be  trusted. 

Note  thaV  in  some  of  the  investigations  reviewed,  the  material  has  been  quenched  to  simulate  a  welded 
condition  or  prestrained  at  high  temperature  to  simulate  creep  damage.  This  usually  leads  to  a  deterioration 
of  the  crack  resistance  and  should  be  considered  when  applying  crack  resistance  values  to  practical  cases. 
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TABLE  1.  CORRELATION  OF  da/dt  WITH  K,  o  .  and  o  , 

*  net  ref 

I  Good  correlation  of  da/dt  with  a  bad  correlation  with  K 

t)  Nicholson  and  Formby  (21):  type  316  stainless  steel,  solution  heat  treated,  SENT  and  CNT  specimens,  740°C. 

2)  Neate  (22):  0.5Cr-0.5Mo-0.25V  steel,  normalized  and  tempered,  SENT  specimens,  565  °C 

To  7  rather  than  o  ). 
ref  •  net 

3)  Neate  (23):  0. 5Cr-0. 5Mo-0. 25V  steel,  normalized  and  tempered,  SENT  and  CT  specimens,  565  °C 

To  7  rather  than  a  ). 
ref  net 

II  Bad  correlation  of  da/dt  with  K,  o  _  not  considered 

- —  net - 

1)  Robson  (24):  0.24  Z  carbon  steel,  normalized,  SEN  and  CT  type  specimens,  450  °C. 

2)  Ellison  and  Walton  (19):  ICr-IMo-V  steel,  SENT  specimens,  565  °C,  correlation  especially  bad  in  bending. 

3)  Haigh  (25):  Cr-Mo-V  cast  steel,  normalized  and  tempered,  WOL  specimens  (ss  CT) ,  550  °C. 

4)  Gooch  (26):  0.5Cr-0.5Mo-0.25  V  steel,  several  quench  rates,  SENB  specimens,  565  °C. 

5)  Gooch  (27):  0.5Cr-0.5Mo-0.25  V  steel,  various  quench  rates  and  tempers,  SENB  and  SENT  specimens. 

6)  Sadananda  and  Shahinian  (28):  Nickel  Alloy  718,  aged,  CT  specimens,  650  °C,  crack  turn.  ling. 

III  Good  correlation  of  da/dt  with  a  K  not  considered 

- net - 

1)  Taira  and  Ohtani  (29):  1Cr-1Mo-0.25  V  steel,  CNRS  and  DENT  specimens,  6000  °C. 

2)  Nicholson  (30):  type  316  stainless  steel,  solution  heat  treated,  DENT  specimens,  600-850  °C. 

IV  -^od  correlation  of  da/dt  with  o  and  with  K 

- net - 

1)  Neate  and  Siverns  (20):  0.5Cr-0.5Mo-0.25V  steel,  quenched,  SENT  specimens,  565  °C. 

2)  Henshall  et  al.  (31):  type  316  stainless  steel,  solution  heat  treated,  tempered  and  untempered, 

DENT  specimens,  636  “C;  1Cr-0.5Mo  steel,  bainite,  SENT  specimens,  565  “C  (°rej  rather  than  o  ). 

3)  Lloyd  (32):  type  316  stainless  steel,  DENT  specimens,  625  °C. 

4)  Taylor  and  Batte  (33):  ICr-Mo-V  steel,  forged,  DENT  specimens,  550  °C,  reinitiation  rather  than 
crack  growth. 

V  Good  correlation  of  da/dt  with  K,  a  not  considered 

- net - 

1)  James  (34):  type  316  stainless  steel,  CT  and  CNT  specimens,  538  °C. 

2)  Thornton  (35):  forged  and  cast  CrMoV  steel,  SENB  specimens,  550  °C. 

3)  Landes  and  Wei  (36):  AISI  4340  type  steel,  CNT  and  constant  K  specimens,  20-140  °C. 

4)  Kenyon  et  al.  (37):  RR  58  aluminium  alloy,  contoured  DCB  specimens  (constant  K) ,  100-200  °C. 

5)  Pilkington  et  al.  (38):  0.5Cr-0.5Mo-0.25V  steel,  vacuum  melted,  SENB  specimens,  550  °C. 

6)  Floreen  (39):  various  nickel  alloys,  CT  specimens,  500-760  °C. 

7)  Kaufman  et  al.  (40):  2219-T851  aluminium  alloy,  CT  specimens,  100-177  °C. 

8)  Worswick  and  Pilkington  (41 ) ;  0.5Cr-0.5Mo-0.25V  steel,  ferritic,  prestrained  at  666  °C.  SENB  SDecimens. 

355^ —  —  — 

9)  Sadananda  and  Shahinian  (28):  type  718  nickel  alloy,  aged,  CT  specimens,  540  °C,  Udimet  700,  aged, 

CT  specimens,  850  °C. 

VI  Bad  correlation  of  da/dt  with  o  K  not  considered 
No  publications. 

VII  Bad  correlation  of  da/dt  with  a  ,  good  correlation  with  K 

• - net  — - 

1)  Siverns  and  Price  (42):  2.25Cr-1Mo  steel,  SENT  specimens,  565  °C. 

2)  Neate  and  Siverns  (20):  2.25Cr-1Mo  steel,  quenched  for  weld  simulation,  SENT  specimens,  565  °C. 

3)  Koterazawa  and  Iwata  (43):  type  304  stainless  steel,  CNRS  and  DENT  specimens,  650  °C. 

4)  Neate  (22):  0. 5Cr-0.5Mo-0. 25V,  quenched  for  weld  simulation,  SENT  specimens,  565  °C. 

6.  CORRELATION  OF  THE  CREEP  CRACK  GROWTH  RATE  WITH  THE 
CRACK  (TIP)  OPENING  DISPLACEMENT  RATE 

For  a  cracked  specimen  as  shown  in  figure  7,  the  displacement  in  the  load  line  v  is  related  to  load  P 
and  compliance  C  as  follows 


Differentiation  with  respect  to  time  gives 


dv  _  dP  .  p  dC  da 
dt  dt  da  •  dt 


For  constant  load 


dv  „  dC  <ia 
dt  da  '  dt 
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For  a  so-called  constant  K  specimen,  ^  is  constant.  It  follows  that  only  for  a  constant  K  specimen  in  the 

elastic  case  at  constant  load  there  will  be  a  straight  proportionality  between  ^  and 
If  these  three  conditions  are  not  fulfilled  a  more  general  relationship 


£-»•(£)' 


may  still  hold. 

It  is  somewhat  surprising  therefore  that  Haigh  (1)  analysing  the  results  of  Robson  (24)  for  a  creeping  CT 
specimen  finds  a  straight  proportionality  nevertheless  (Fig.  14).  As  shown  in  figures  15  and  16  the  time 
dependence  of  COD  and  crack  length  a  can  be  quite  different.  A  relation  as  above  will  have  only  a  limited 
validity. 

For  a  materiel  behaviour  as  depicted  in  figure  16  it  may  be  stated  that  for  the  onset  of  crack  growth  v 
must  first  reach  a  critical  value  v^,  and  that  the  relationship  will  than  hold  for  da/dt  and 

4-  (v  -  v  )  (6.5) 

dt  c 

If  the  relationship  does  not  hold  this  may  show  up  as  a  non-el iminated  dependence  on  load  P  Pig.  17). 
Radhakrishnan  and  McEvily  (9,  44)  found  that  in  that  case  the  data  points  could  be  made  to  fall  on  one 
•  •  •  del 

straight  line  by  plotting  ^  versus  the  parameter 


mr 


They  claim  a  theoretical  background  for  the  relationship 
da  ,  ,  f  dv  /PcAaTt3 


o[s(?)7 


where  a  and  0  are  functions  of  temperature  T  (45).  Sometimes,  as  for  Ti-6Al-2Sn-4Zr-2Mo-0 . 1  Si,  a  will  be 
very  small  and  even  practically  zero  (46). 

While  v  and  dv/dt  are  directly  related  to  the  overall  energy  dissipation,  a  more  mechanism  oriented 
approach  makes  use  of  the  crack  tip  opening  displacement  6.  The  trouble  is  that  6  is  much  more  difficult  to 
determine  than  v. 

Nicholson  (30,  47)  and  Henshall  et  al.  (31)  determined  6  on  metal lographic  specimens  made  after  the 
cracking  test.  They  correlated  &  with  v  and  then  determined  ^  from  the  time  record  of  v. 

The  more  usual  method  is  to  assume  that  there  is  a  plastic  hinge  at  a  certain  spot  in  the  unbroken  ligament 
b  and  that  6  will  be  a  certain  fraction  of  v. 

Haigh  (25)  found  the  relationship 


da  _  „  / d6\^ 

dt  =  1  ‘  \dt,/ 


Other  investigators  such  as  Pilkington  et  al.  (38),  Nicholson  and  Formby  (21),  Gooch  et  al.  (26,  48)  and 
Henshall  et  al.  (31)  have  found  relations  of  the  type 

(a  -  aQ)  =  1  (6  -  60)  (6.9) 

where  u  is  similar  to  the  crack  tip  opening  angle. 

The  question  is,  given  a  relationship  between  and  4r  or  4r»  what  the  practical  use  of  it  is.  There 
_  ,  at  at  at 

are  two  needs:  ^ 

-  for  a  high  temperature  installation  one  would  like  to  determine  ^  from  a  more  easy  to  measure  quantity  or 

-  one  would  like  to  predict  from  a  more  easy  to  calculate  quantity. 

The  conclusion  is  that  the  relationship  will  be  of  little  practical  use  because  for  a  structure  and 

da  ^ 

are  neither  more  easy  to  measure  nor  more  easy  to  calculate  than  -j—  itself. 

7.  CORRELATION  OF  THE  CREEP  CRACK  GROWTH  RATE  WITH  THE 
J- INTEGRAL  OR  THE  RELATED  PARAMETER  J* 


In  section  3.3  the  J-integral  has  been  defined  already  and  a  method  was  indicated  to  determine  J 
experimentally  (Fig.  7). 

With  regard  to  creep  cracking  Landes  and  Begley  (15)  have  applied  a  parameter  fhat  strongly  resembles  the 
J-integral.  They  have  indicated  it  as  C*  whereas  other  authors  use  the  symbol  J.  To  avoid  confusion  only 
the  symbol  J*  will  be  used  in  the  present  paper.  The  integral  in  question  is  defined  as  (see  Fig.  18) 

J*  -  fa  |w*dy  -  Tj  dsj  (7.1) 

where  c 

mn 

W*  -  f  o.kdrjk  (7.2) 


It  appears  that  J*  can  be  obtained  from  J  by  substituting  strain  rates  for  strains  and  displacement  rates 
for  displacements.  It  is  clear  that  in  general  J*  ¥  dJ/dt  and  that  hence  the  indication  J  is  misleading. 
Like  J,  J*  can  be  calculated  by  finite  element  analysis  or  from  analytical  stress-strain  relations. 

For  J  an  empirical  determination  was  possible  on  account  of  the  relation 
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_  _  dU  _ 

- d  I 

da 

da  \ 

Likewise  J*  can  be  determined  from 


H'~} 

determined  fr 

-U'H 


Landes  and  Begley  determined  J*  from  cracking  tests  at  constant  displacement  rate  v.  They  processed  their 
test  data  as  shown  in  figure  19.  Diagrams  are  produced  in  the  following  sequence 

-  Load  and  crack  length  as  a  function  of  time 

-  crack  growth  rate  as  a  function  of  crack  length  at  constant  displacement  rate 

-  load  as  a  function  of  crack  length  at  constant  displacement  rate 

-  load  as  a  function  of  displacement  rate  at  constant  crack  length 

-  energy  rate  U*  as  a  function  of  crack  length  at  constant  displacement  rate 

-  energy  rate  derivative  J*  as  a  function  of  displacement  rate  (at  constant  crack  length) 

-  crack  growth  rate  (at  a  certain  displacement  rate  and  crack  length)  as  a  function  of  J*. 

Sadananda  and  Shahinian  (49)  have  developed  a  method  to  determine  both  J  and  J*  from  the  results  of  creep 
cracking  tests  at  a  series  of  constant  loads  P.. 

Their  schemes  are  shown  in  figure  20  which  is  self-explanatory. 

Many  authors  have  used  approximate  expressions  from  which  J  and  J*  are  more  easy  to  determine.  These  will 
be  discussed  later. 

The  question  posed  now  is  whether  there  is  an  unambiguous  relation  between  da/dt  and  J  or  J*,  i.e.  a 
relation  that  will  hold  even  for  different  geometries  and  sizes  of  specimens,  different  load  levels  and 
initial  crack  lengths  etc.  Only  in  that  case  will  one  be  able  in  principle  to  accurately  predict  da/dt  from 
J  or  J*. 

Landes  and  Begley  (15)  testing  CCP  (=CNT)  specimens  of  the  FeNiCrMo  alloy  Discalloy  at  647  °C  found  a  single 
relation  between  da/dt  and  J*  if  the  same  specimen  was  tested  either  at  a  single  or  at  a  series  of  dis¬ 
placement  rates  v.  If,  however,  they  tested  CT  specimens  rather  than  CCP  they  found  the  curve  to  be  offset 
relative  to  the  former  (Fig.  21).  The  correlation  was  better,  however,  than  with  K  or  nominal  stress 
(■  reference  stress). 

Nikbin,  Webster  and  Furner  (50)  tested  the  aluminium  alloy  RR58  at  150  °C  and  steel  type  0.5Cr-0.5Mo-0.25V 
at  565  °C.  They  used  contoured  DCB  specimens  which  give  a  constant  K  at  constant  load  P.  They  found  that 
da/dt  was  not  constant.  If,  however,  they  adjusted  P  repeatedly  to  arrive  at  a  constant  J*  they  found  that 
da/dt  was  constant  for  a  large  interval  of  crack  lengths,  i.e.  excluding  small  primary  and  a  tertiary 
regions. 

Harper  and  Ellison  (51)  tested  ICrMoV  steel  at  565  °C.  They  used  SENB  and  CT  specimens.  The  correlation 
of  da/dt  with  J*  was  quite  good  at  high  values  of  these  quantities.  Curves  for  constant  bending  moment 
were  fanning  out  however  at  low  values,  giving  lower  da/dt  values  at  higher  applied  moments.  The  authors 
ascribed  it  to  the  effect  of  primary  creep. 

Sadananda  and  Shahinian  (49)  tested  nickel  alloy  718  at  538,  649  and  760  “C  using  CT  specimens.  They 
found  that  for  a  series  of  applied  loads  P,  the  crack  growth  rate  correlated  rather  well  with  both  J  and 
J*  and  rather  badly  with  both  stress  intensity  factor  K  and  nominal  or  reference  stress  o  . 

Pilkington  and  Worswick  (52)  continuing  work  reported  in  (41)  tested  0.5Cr-0.5Mo-0.25V  atn649  °C  using 
SENB  specimens.  The  material  had  been  prestrained  plastically  0.26  resp.  0.58  X  at  666  °C  to  introduce 
cavitation  damage.  For  different  applied  loads  they  found  that  da/dt  correlated  rather  well  with  d6/dt, 
and  with  quantities  J*(v)  and  J*(6)  defined  as 


.....  P  dv  ,  P  dd 

J*(V)  "  ta  •  di  3nd  J*(6)  Ba  *  dt 


(7.5,  7.6) 


Earlier  they  had  reported  good  correlation  with  log  K  but  had  found  distinct  effects  of  primary  creep 
and  prestraining.  This  no  longer  showed  up  in  their  newer  plots.  This  will  be  because  cavitation  not  only 

increases  ^  but  also  and  and  the  same  may  be  true  for  primary  creep. 

Branco  and  Radon  (53)  have  continued  the  work  of  Kenyon  et  al.  (37)  trying  to  correlate  da/dt  not  with  Y 
but  with  J.  They  did  find  a  correlation  but  did  not  report  whether  this  was  any  better  than  with  K, 

a  .or  J*.  They  introduce  the  relation 
net 

(7-7) 

where  C  and  J  are  strongly  dependent  on  temperature  T.  Since  the  relations  with  T  involved  are  not 
monotonic  the  validity  of  the  equation  is  doubtful. 

Saxena  (54)  tested  annealed  304  type  stainless  steel  at  temperatures  between  538  and  705  °C.  He 
applied  CCT  and  CT  specimens  using  constant  displacement  rates  v  that  were  sometimes  varied  in  a  stepwise 
manner.  The  correlation  of  da/dt  with  J*  was  rather  good  but  deviations  occurred  for  very  long  cracks  in 
the  CCT  specimens. 

In  Japan  very  comprehensive  investigations  have  been  made  to  determine  the  validity  of  J*  as  a 
governing  parameter  for  creep  crack  growth. 

Koterazawa  and  Mori  (55,  56  and  Fig.  22)  tested  type  304  stainless  steel  at  650  °C  in  air.  They  used  three 
different  double  edge  notch  specimens  and  one  central  notch  specimen. 

The  material  came  from  three  heats  (A,  B  and  C)  which  was  accounted  for  in  the  specimen  designation, 
e.g.  DEN  2C.  They  approximated  J*  as 

J*  ■  °net  *  £  <7‘8> 

The  crack  growth  rate  appeared  to  correlate  rather  well  with  J*  (Fig.  22).  The  correlation  with  K  was  quite 
bad,  the  one  with  o  was  much  better,  and  better  still  if  a  was  divided  by  w  to  eliminate  a  size  effect 
(Fig.  23).  net 

Ohji  (4)  and  Ohji  et  al.  (57)  tested  type  304  stainless  steel  at  650  6C  using  a  large  number  of  specimens 
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of  various  sizes  and  shapes.  Figures  24  and  25  show  that  the  correlation  of  da/dt  with  J*  was  superior 
to  those  with  K  and  o  .  This  is  quite  obvious  also  for  some  added  results  for  CT  specimens. 

In  (57)  results  are  p?lsented  also  for  circumferentially  notched  round  specimens  at  three  diameters. 

Results  have  been  added  of  tests  on  single  and  double  edge  notch  specimens  as  well  as  central  notch 

specimens . 

In  a  diagram  of  da/dt  versus  J*  all  results  fall  in  a  relatively  narrow  scatter  band. 

Taira,  Ohtani  and  Kitamura  (58)  have  done  comprehensive  tests  on  0.16  %  carbon  steel  at  400  and  500  °C  in 
air,  at  400  bC  in  vacuum,  on  type  304  stainless  steel  at  650  QC  in  air  and  in  vacuum,  on  type  316  stainless 
steel  at  600  and  650  °C  in  air. 

They  used  three  types  of  specimens,  each  in  three  different  sizes  (Large,  Medium  and  Small)  i.e. 

-  hollow  cylindrical  specimens  with  a  local  through  crack:  CNC-L,  CNC-M,  CNC-S. 

-  solid  cylindrical  specimens  with  a  circumferential  notch:  RNB-L,  RNB-M,  RNB-S. 

-  flat  plate  specimens  with  a  central  notch:  CNP-L,  CNP-M,  CNP-S. 

J*  was  again  approximated  as  onet  • 


As  figure  26  shows,  the  J*  integral  approach  can  cope  quite  well  and  much  better  than  the  c 
with  differences  in  stress  level  and  specimen  size  and  shape. 


approach 


An  investigation  similar  to  those  conducted  in  Japan  was  carried  out  in  Great  Britain  by  Christian 
et  al.  (59).  These  authors  tested  0.5CrMoV  steel  and  ICrMoV  steel  at  565  °C  using  CT,  SENB  and  DCB  type 
specimens.  For  these  different  geometries  they  tried  to  correlate  the  crack  growth  rate  with  K,  cr  and  J*. 
Correlations  with  J*  were  more  satisfactory  throughout.  The  diagrams  they  have  produced  closely  resemble 
those  by  the  Japanese  authors  shown  in  figures  24  and  25,  meaning  that  in  the  diagram  of  da/dt  versus  J* 
there  was  still  quite  a  bit  of  spread  although  it  was  less  than  with  the  other  parameters. 

Seemingly  alarming  results  in  this  context  have  been  presented  by  Donath  et  al.  (60).  These  authors  tested 
the  nickel-base  superalloy  IN  100  at  732  “C  using  two  different  specimen  types.  The  first  one  was  the 
comnon  CT  specimen  but  the  second  a  closed  ring  tested  in  tension  and  cracked  from  the  inner  side  at 
locations  90°  away  from  the  points  of  load  application.  A  range  of  thicknesses  was  used  for  both  specimen 
types.  They  evaluated  the  stress  intensity  factor  K,  net  section  stress  a  and  the  contour-integral  J* 
as  correlating  parameters  describing  crack  growth.  A  special  feature  of  their  investigation  was  that 
cracking  occurred  at  increasing  K  in  the  case  of  the  CT  specimens  and  at  practically  constant  K  in  the  case 
of  the  ring  specimens.  They  applied  a  scheme  like  that  in  figure  20  to  determine  J*  from  the  test  results. 
They  encountered  a  problem  already  encountered  by  other  authors  (e.g.  Worswick  and  Pilkin^ton,  41)  in  that, 
due  to  primary  creep  effects,  crack  races  were  high  initially,  gradually  decreased  to  a  minimum  and  then 
increased  again  at  increasing  crack  length  or  remained  constant.  This  means  that  in  the  da/dt  versus  K 
diagram  a  V-shaped  or  root-sign-shaped  curve  (Fig.  28  top)  is  obtained  rather  than  the  sigmoidal  one 
sketched  in  figure  9. 

Barring  the  initial  downward  branches  Donath  et  al.  found  good  correlations  in  spite  of  thickness 
variations  of  da/dt  with  all  three  parameters  K,  anet  and  J*  (Fig.  27).  The  situation  was  different  for  the 
ring-shaped  specimen,  however. 

In  the  da/dt  versus  K  plot  they  obtained  a  series  of  parallel  vertical  lines  corresponding  to  the  various 
initial  K-values  (Fig.  28).  The  thickness  effect  was  negligible.  Since  during  each  test  K  was  practically 
constant,  each  line  depicted  nothing  but  the  decreasing  crack  growth  rates  due  to  primary  creep.  This 
situation  really  is  not  alarming  since  the  bottom  ends  of  the  lines  lay  fairly  close  to  the  da/dt  versus  K 
curve  derived  for  CT  specimens.  They  considered  the  onet  approach  useless  for  the  ring  specimen,  but  then 
it  ought  to  be  possible  to  define  a  reference  stress  o  f  even  for  this  geometry.  They  claim  the  J* 
approach  to  be  useless  too  for  the  ring  specimen  because  ofrthe  decreasing  crack  growth  rate  at  increasing 
crack  length.  Apparently  they  have  not  realized  that  the  J*  approach  should  be  restricted  to  secondary 
creep  and  crack  growth  rates  and  excludes  the  primary  ones. 

All  in  all  the  impression  is  gained  that  the  J*  integral  holds  better  promise  as  a  governing  parameter  than 
K  or  onet  or  °re£t  although  the  situation  is  still  not  ideal. 

One  deterrent  to  working  with  J*  is  the  laborious  way  in  which  it  has  to  be  determined  if  one  strives  for 
greater  accuracy  than  with  the  approximations  used  e.g.  by  the  Japanese  authors.  Fortunately  several 
alternative  methods  are  available  for  the  quantitative  determination  of  both  J  and  J*.  These  will  be 
discussed  in  the  next  sections. 

8.  ALTERNATIVE  EMPIRICAL  DETERMINATION  OF  J  AND  J* 


In  figures  7  and  20  methods  have  been  shown  to  determine  J  from  a  series  of  load-displacement  diagrams. 
Rice  and  other  authors  have  developed  methods  to  determine  J  from  one  single  load-displacement  diagram, 
i.e.  from  one  test. 

In  principle  the  same  methods  can  be  used  for  J*  working  with  a  diagram  of  load  versus  displacement  rate, 
but  it  would  require  more  than  one  test  to  obtain  a  load-displacement  rate  diagram. 

J  has  been  related  to  the  energy  dissipation  as 

j.-iM 

where  B  is  specimen  thickness  and  U  is  the  area  below  the  load  displacement  curve  (Fig.  29,  A  and  B) . 

Now 

v  P 

U  -  /  Pdv  •  Pv  -  /  vdP  -  Pv  -  U  (8.2) 

0  0 

It  follows  that: 


i  j7.  *l) 

B  0  V  9a'v 
i  *(—)  dP 

»  n  Wp 


J-i0'^pdP 

It  can  be  shown  that  the  latter  two  integrals  can  be  build  up  from  elementary  areas  of  the  P-v  diagram 
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such  as  shown  in  figure  29C. 

According  to  Rice  (61)  it  has  advantages  to  split  up  v  in  two  components. 
Either 


v  +  v 
noc  r  c  r 


v  ,  +  v  , 
el  pi 


In  the  first  case  v  has  been  split  up  in  a  part  that  would  occur  in  the  absence  of  a  crack  and  in  a  part 
that  is  contributed  by  the  crack. 

In  the  second  case  v  has  been  split  up  in  elastic  and  plastic  contributions. 

For  a  flat  double  edge  crack  specimen,  where  the  unbroken  ligament  measures  b,  loaded  by  a  force  per  unit 

thickness  P,  it  can  be  shown  that  v  ,  should  have  the  form 

Pi 

vpl  *  b  “(I)  (8'7) 

where  H  is  an  unknown  function  of  (P/b) . 

It  can  be  derived,  however,  that 
p/  3v  \ 

J  -  Jel  +  /(-  -&■)«  (8'8> 


It  follows  that 


P  ■  f  "’({)  ■<{)•¥  K^l  -  v} 


Moreover,  J  is  equal  to  the  elastic  strain  energy  release  rate  G  -  K  /E.  Partial  integration  leads  to 


My 


Now  the  quantities  between  parentheses  can  be  identified  as  elementary  areas  in  figure  29C. 

v  ,  is  usually  neglected, 

pl.nocr  1  a 

It  is  seen  that 


Pv  ,  »  Pv  .  »  U  ,  +  U  , 

pi  pl.cr  pl.cr  pl.cr 


and  that 


r  pl 

J  Pdv 


-  U  .  -  Pv  ,  -  U  , 

pl  pl.cr  pl.cr  pl.cr 


(8.12) 


J-G+iju.  -  U  ,  1 

b  \  pl.cr  pl.cr J 


(8.13) 


The  same  results  holds  for  a  plate  with  a  central  crack. 

Rice  et  al  have  made  similar  analyses,  leading  to  similar  but  different  expressions  for  a  circumferentially 
notched  round  specimen,  for  a  cracked  three-point  bend  specimen  and  for  a  compact  tension  specimen. 

In  the  last  case  there  is  found 


J  -  |  J  Pdv  }  w  £  U 

b  J  cr  b  \  tot  nocrj  b  t 


These  and  other  analyses  have  been  reviewed  by  Keller  and  Munz  (62). 

Kanazawa  et  al.  have  attempted  to  derive  more  accurate  expressions  for  the  three-point  bend  specimen  and 
the  compact  tension  specimen. 

Merkle  and  Corten  have  derived  a  more  complicated  and  supposedly  more  accurate  expression  again  for  the 
compact  tension  specimen. 

In  a  more  recent  publication  Clarke  and  Landes  (63)  have  simplified  the  Merkle  and  Corten  analysis  for  the 
compact  tension  specimen. 

They  arrived  at  the  expression 

<*•'» 

Here  a  is  a  rather  complex  expression  in  crack  length  a  and  ligament  b.  However,  the  authors  present  an 
approximation  obtained  by  curve  fitting 


1  ♦  0.261 


('  -  i) 


(8.16) 


Rice  et  al.  (61)  have  derived  also  equations  with  which  J  can  be  calculated  from  one  single  point  on  the 
P-v  curve.  The  accuracy  is  doubtful  however,  and  the  result  may  serve  only  as  a  quick  estimate. 

9.  ANALYTICAL  DETERMINATION  OF  J  AND  J*  USING 
NON-LINEAR  STRESS-STRAIN  RELATIONS 


Rice  and  Hutchinson  and  several  of  their  co-workers  have  calculated  stresses  and  strains  around  the 


tip  of  a  stationary  crack  on  the  basis  of  a  quasi-plastic,  non-linear  elastic  relation  between  stress  and 
strain 

\N  _  /_ 


%  *  ft/  °r  '  “  (%) 


and  accounting  for  tri-axiality  by  means  of  the  Von  Mises  effective  stress 

o  2  *  4  S.'.  S.  .  (9.2 

e  2  ij  ij 

where  S. .  are  the  deviatoric  stresses. 

J-J 

By  working  with  polar  coordinates  around  the  crack  tip,  the  J-integral  can  be  calculated  in  a  straight¬ 
forward  manner  by  integrating  along  a  circular  path. 

It  should  be  noted  that  in  these  analyses  so-called  deformation  theory  of  plasticity  has  been  used, 
which  makes  the  results  valid  only  for  monotonically  increasing  stresses. 

By  replacing  the  above  stress-strain  relationship  by  the  so-called  Norton  law 


the  J*-integral  can  be  calculated  in  an  identical  way  by  making  use  of  the  elastic  analogue,  although  the 
equations  have  been  derived  for  a  non-growing  crack. 

The  publications  with  which  it  all  started  are  Rice  and  Rosengren  (64)  and  Hutchinson  (65,  66). 

Hilton  and  Hutchinson  (67)  present  the  results  by  introducing  plastic  stress  and  strain  intensity  factors 


o.  .  -  K  .  r 
ij  o 


.  a«>^j(0) 


e.  .  “  K 
ij  e 


1 +m  / . . 

r  .  e»>^j(0) 


where  and  are  directly  related  to  the  elastic  stress  intensity  factor 
2 

K  -  c  (K  .//if) 1+m 
o  m  el 


c»  ‘  (t)1 


Im  is  a  rather  complicated  integral  tabulated  by  Hutchinson  (66). 

The  authors  have  applied  the  analysis  to  a  variety  of  specimen  shapes  and  for  the  cases  of  plane  strain  as 
well  as  plane  stress.. 

Goldman  and  Hutchinson  (68)  write 


J  -  aone-K  K  I  (9.9 

0  0  o  e  m 

The  analysis  is  applied  to  a  flat  strip  of  width  2w  and  a  central  crack  2a,  loaded  in  plane  strain  by  a 
stress  0“  at  infinity.  They  also  calculate  the  COD  A  in  the  axis  of  symmetry.  The  results  are  presented 

as  functions  of  ~  and  m  and  they  are  normalized  with  respect  to  the  parameters  from  the  stress-strain 

relation,  crack  length  and  applied  stress  o°°  or  strain  c”. 

The  equations  are  as  follows 

A  *  AN(t-m)(2EOa  ■  ANfe*mX2eOaa  T-  (r  (9- 1 '  > 

J  is  normalized  also  with  respect  to  A 


/a  NT  - 1 /a  | 

f  w  'j 

1  /m 
I  | 

f  4  ^ 

t 

1  L 

\w-7T  °oeoa| 

U  -  *) 

1  1 

V2W 

1  / 

Goldman  and  Hutchinson  have  calculated  also  the  load  line  displacement  v,  and  have  split  it  up  as 

v  •  v  ♦  v  where  v  -  2hE»  and  2h  is  the  strip  length, 
cr  nocr  nocr  r  0 

The  results  are  presented  as 


vcr  ‘  vcrN  *  A 


The  quantities  and  A^  appear  to  depend  strongly  on  Therefore  they  are  normalized  further  with  respect 
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to  a  and  w. 

The  following  forms  have  been  chosen 


'  ■  fji 

A  *  ANaw(^“1){2£0a(^l)  *  0  -  ANaw(w’m){2£Oaa(w_^_a)  T  (f 


(9.14) 

(9.15) 


The  numerical  results  are  presented  in  table  2. 

Shih  and  Hutchinson  (69)  have  continued  the  work.  They  consider  a  flat  strip  with  central  crack  in  plane 
stress.  They  normalize  their  results  with  respect  to  load  P  per  unit  thickness  and  a  reference  value  P^ 


J  -  g,(^m){aa0e0a(l  -  ^)(|-)  }  (9.16) 

A  _  g2(f«m){a£0a(y  }  (9-,7) 

vcr  ‘  *3(>){a£0a(y  }  (9’18) 

Note  that  P  -  2wo“  and  that  ■  — 7^ - 7  (9.19,  9.20) 

pn  °n(“  "  a' 


Numerical  values  are  given  in  table  3. 

The  results  in  their  present  form  do  not  allow  interpolation  for  m  >  10  and  to  the  limit  values  “  ”  0  and 

^  -  1.  An  adapted  way  of  normalizing  is  suggested  by  Amazigo's  exact  solution  for  the  anti-plane  shear 
case.  This  leads  to  the  equation 


+  t(n 

-  1) 

a 

W 

-1 

rrrij 

^  m) 

+  81 

(b*  J 

J 

m 

(9.21) 


The  result  is  that  G(— ,m)  varies  far  less  strongly  with  —  and  m  than  g.(— ,m)  and  thereby  allows  extra¬ 
polation.  w  w 

Shih  and  Hutchinson  state  that  the  calculations  presented  so  far  are  good  approximations  for  the  case  of 
v,large  scale  yielding"  where  the  contributions  from  elasticity  can  be  be  neglected. 

For  the  elastic-plastic  case  they  propose  to  add  the  linear  elastic  to  the  non-linear  elastic  solution. 
The  linear  elastic  solution  is  of  course  a  special  case  of  the  non-linear  elastic  one. 

Symbolically 


J  »  J(a  -  ae,m  -  1)  +  J(a,m)  (9.22) 

Note  that  they  apply  the  Irwin  correction  in  the  elastic  case,  i.e.  ag  ■  a  +  r^. 

In  the  present  notation 


rY_1  m-1  w-a  /a  A/P  \ 

a  2it  '  m  +  1  a  '  81\w’  /\pq/ 

(9.23) 

This  leads  to  the  following  equations 

2  m+1 

w(«T1)’*,8,(W,,Xpo)  +a8’tmX0 

(9.24) 

(9.25) 

v  *  ^ 

(9.26) 

Here: 

a  *  a  +  r  for  P  £  P. 
e  Y  u 

a  »  ( a  )  for  P  >  P_ 

r  -  1  m  f  ^ 

Y  2tt  m  ♦  1  \Oq/ 

(9.27) 

(9.28) 

(9.29) 

a 

,  e  w-a 

*  ■  —  . - 

a  w-a 

(9.30) 

e 


Results  of  calculations  using  these  equations  are  reported  to  agree  well  with  results  of  FEM  calculations. 
Shih  and  Hutchinson  (69)  have  considered  also  the  single  edge  crack  strip  loaded  in  bending  in  the  plane 
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strain  condition.  The  width  of  the  uncracked  ligament  is  b  •  w  -  a. 

For  reasons  of  brevity  only  the  equations  for  the  elastic-plastic  case  are  presented  here. 


^I-*2I£h2(^,)^  +  ah2(>)(ym 
^  •  ♦SGr’1)  t0  * 


(9.31) 


(9.32) 


(9.33) 


Note  that  8^^  is  the  contribution  of  the  crack  to  the  relative  rotation  of  the  ends  of  the  strip. 
Some  results  of  a  FEM  analysis  for  —  •  0.5  are  presented  below. 


m  “  1 

m  “  2 

m  ■  3 

m  ■  5 

m  -  7 

m  -  10 

1,104 

0,957 

0,851 

0,717 

0,653 

0,551 

5,129 

3,640 

2,947 

2,255 

1 ,953 

1,606 

2,749. 

2,359 

2,032 

1,590 

1,373 

1,121 

Also  in  this  case  there  was  good  agreement  between  results  of  analytical  and  of  FEM  calculations. 


stress  and  in  plane  strain. 

They  present  the  equations 

h1(in){%c0b(|i-)"  '} 

(9.34) 

A  -  h2(>){6i0a(f-)  } 

(9.35) 

*cr  ■  ^)M^} 

(9.36) 

Here: 

Pg  *  1 ,455  nbOg  (plane  strain) 

(9.37) 

Pq  *  1,072  nbOjj  (plane  stress) 

(9.38) 

(9.39) 

Results  of  their  calculations  are  presented  in  tables  4  and  5.  In  another  publication  (71)  the  authors  give 
the  same  tables  but  then  relating  to  time  ^independent  plastic  deformation.  In  the  latter  publication  they 
also  treat  the  elastic-plastic  case  in  the  way  discussed  before. 

An  interesting  question  is  whether  values  for  J  and  J*  calculated  from  the  above  analyses  agree  with 
values  determined  from  P-v  diagrams  or  P-v  diagrams. 

Saxena  (34)  has  compared  the  results  from  both  techniques  using  creep  law  parameters  reported  in  the 
literature.  He  found  discrepancies  between  the  derived  J*  values  to  range  from  8  to  40  Z,  which  is  not  bad 
at  all,  considering  the  many  possible  sources  of  error. 

10.  VARIOUS  OTHER  MAYS  TO  CALCULATE  J  AND  J* 

Nikbin,  Webster  and  Turner  (30)  have  calculated  J  and  J*  for  a  Contoured  Double  Cantilever  Beam 
spec imen . 

They  assume  the  non-linear  stress-strain  relation 


Webster  has  shown  that  for  a  material  that  obeys  this  law  J  can  be  calculated  as 


.  P  dv 

B  (m  +  1)  '  da 
n 


(10.1) 


(10.2) 


The  quantity  ^  is  calculated  for  a  cantilever  beam  of  length  a  and  the  final  result  is 

.  2<ae0°o‘B)  f 2m  ♦  llm  (aP)m*' 

Bn<n  ♦  1)  •  llSB- /  •  (h/2)2m+1 


(10.3) 


Here  B  and  B  are  the  nominal  thickness  and  the  local  thickness  in  the  side  groow  and  h  is  the  local  beam 
height.  n 

For  a  material  that  obeys  Norton's  law 
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i  -  e(— )" 

s0  '°0' 


there  is  found  similarly 


_  2(§£0°0  ^  J (2n  +  1)l"  (aP)r 
“B7H-rTTl  2nB  J  ,hm2r 


These  equations  will  not  be  very  accurate.  Similar  calculations  have  been  made  for  stress  corrosion 
specimens  strained  elastically  and  important  errors  were  found  to  result  from  ignoring  the  deformation  by 
shear  and  the  lack  of  rigidity  at  the  "clamped"  cross  section  at  the  crack  tip. 

Branco  and  Radon  (53)  have  refined  the  above  analysis  by  considering  elastic  and  creep  deformation  in  both 
bending  and  shear.  They  also  account  for  the  crack  tip  influence  by  adding  an  extra  elastic  deformation  in 
bending  calculated  by  assuming  that  the  beam  is  clamped  at  a  distance  i  b  from  the  crack  tip  inside  the 
ligament  b. 

Nikbin  et  al.  (72)  give  an  approximate  solution  for  J*  starting  with  the  equation 


B  (n  +  1)  da 
n 


They  assume  that 


v  *  j  .  an  .  g(P) 


so  that 


and  thus 


n 

Harder  and  Ellison  (51)  calculate  J*  using  a  method  based  on  the  reference  stress  concept.  They  assume  that 
during  creep  the  stress  distribution  is  the  same  as  at  limit  load.  This  enables  then  to  calculate  the 
quantity  U*  (see  section  7),  which  after  dif ferentation  with  respect  to  a  gives 


i*  .  ~n  El  fl  <ju  1 

n  ♦  1  ‘  Bw  '  [p  d(a/w) J 


(10.10) 


For  the  usual  large  values  of  n  this  approximates  to 


_  Pv  h  du  1 
Bw  [u  ’  d(a/w) J 


(10.11) 


For  bend  specimens  there  is  found  similarly 


T*  ,  M0  fi  du _ ] 

Bw  [u  ’  d(a/w) J 


(10.12) 


Here  u  is  the  ratio  of  the  yield  load  for  a  cracked  specimen  to  that  for  an  uncracked  specimen.  Values 
of  nas  a  function  of  (a/w)  have  been  compiled  for  a  variety  of  specimens  by  CEGB  in  the  U.K. 

Koterazawa  and  Mori  (55)  start  with  an  equation  for  J*  analgous  to  one  for  J  used  by  Rice  (section  8) 


f 


(10.13) 


In  further  analogy  they  write  for  a  double  edge  crack  specimen 


v  -  bf (P/Bb) 


(10.14) 


and  on  the  basis  of  Norton's  law: 


v  -  Cb(P/Bb) 


(10.15) 


where  C  is  a  constant. 


From  the  first  equation  they  find 


J* 


1  P  /8v\ 
n  +  1  B  \3a/. 


(10.  16) 


P«const. 


and  finally 

n  -  1  Pv 


J* 


n  +  1  Bb 
For  large  n  this  reduces  to 


(10.17) 


r*  PV 

J*  '  Bb 


(10.18) 


This  equation  holds  also  for  a  centre  crack  flat  specimen. 

Similar  equations  have  been  derived,  for  a  circumferentially  notched  round  specimen 


J*  «  Pv/nr 

(where  2r  is  the  diameter  of  the  ligament)  and  for  a  compact  tension  specimen 
2Pv 


(10.19) 


J*  « 


Bb 


(10.20) 


Taria  et  al  (58)  use  similar  equations  developed  by  Ohji. 
For  the  centre  crack  flat  specimen  they  write 


n  -  1  • 

j*  «  — — — •  a  A 


n  +  1  net 

and  for  the  circumferentially  notched  round  one 


(10.21) 


J*  -  ■  1  0  _  A 


2n  +  1  net 


(10.22) 


Here  A  is  the  COD  rate  at  the  axis  of  symmetry  in  the  first  case,  and  the  COD  rate  at  the  outer  surface 
in  the  second  case. 


Mu8icco  et  al.  (73)  have  compared  various  different  methods  to  determine  J*.  These  include 
Hutchinson  approach  leading  to  an  equation  of  the  form 

n+1 


the 


The  Landes  and  Begley  graphical  method  based  on 


(10.23) 


J* 


dU* 


da 


(10.24) 


the  Ohji  et  al.  simplified  equation  (see  also  eq.  10.17) 
P 


J*  ‘  «(n)  v 

the  Harper  and  Ellison  formula  (see  also  eq.  10.10) 


(10.25) 


J*  -  h(n,y)  |-  v 


(10.26) 


and  an  equation  derived  by  Musicco  and  Bernasconi 


J*  -  g(n)  .  oQ  .  S  (10.27) 

They  applied  these  to  CT  specimens  of  304  type  stainless  steel  tested  at  550  °C  in  argon.  They  observed  a 
fair  agreement  for  the  various  formulae  and  with  test  results  reported  by  Saxena . 

11.  THEORETICAL  PREDICTIONS  OF  THE  CRACK  PROPAGATION  RATE 

Various  authors  have  derived  equations  relating  the  crack  propagation  rate  to  a  parameter  such  as  K, 


G,  o  or  J*. 

Theseanalyses  do  not  give  proof  that  Che  parameter  in  question  really  determines  the  crack  propagation 
rate.  Often  assumptions  are  made  that  lead  to  the  result  envisaged. 

For  the  sake  of  completeness  the  formulae  will  be  presented.  For  their  derivation  the  reader  is  referred 
to  the  original  publications. 

Landes  and  Wei  (36)  start  with  Krafft's  ligament  model  and  apply  equations  derived  by  Rice,  Rosengren  and 
Hutchinson  (section  9).  They  arrive  at  the  formula 


2  Nil  2N_ 


Here  K  is  the  critical  stress  intensity  factor  and  c  .  is  the  secondary  creep  rate  in  a  point  at  a 
distance  d  from  the  crack  tip.  m 


f  K  1 

lo' 


75  Ncy) 


Taira  and  Ohtani  (29)  have  made  a  phenomenological  analysis  involving  grain  boundary  creep  rate,  crack 
nucleation  rate  and  fracture  strain,  and  a  typical  length  p*  measured  from  the  crack  tip. 

The  result  is 

if  "  bgbc1c2(pA)"^  0  8b<lx  (,1-3) 

Here  b  .  ,  c.,  c„  and  a  ,  are  constants, 
gb’  1’  2  gb 

Barnby  (74,  75)  applies  Hoff's  elastic  analogue  and  Norton's  law  to  calculate  stresses  near  the 
crack  tip.  He  assumes  the  crack  propagation  rate  to  be  proportional  with  the  creep  rate  at  a  characteristic 
distance  frpm  the  crack  tip.  He  arrives  at  the  formulae 


e'('2m  ~  1 

0\  2m 


K"  1 

2*)°*5J 


(w  -  a) 


Here  is  an  empirical  constant  and  K"  is  a  stress  intensity  factor  under  creep  conditions. 

Barnby  and  Nicholson  (76)  use  Neuber's  rule  in  addition  to  Hoff's  elastic  analogue  and  Norton's  law 
to  calculate  stresses.  They  take  the  crack  propagation  rate  to  be  proportional  to  the  deformation  rate  of 
an  element  of  length  L  orthogonal  to  the  crack  and  at  a  distance  d  from  the  critical  tip.  This  leads  to 

the  result  „ 

2n 


da  .  J  OP  f 

dt  i(EB2wd  )°'5J 


Ohji  (4)  and  Kubo  et  al.  (77)  consider  a  creep  crack  loaded  in  mode  III  (anti-plane  shear).  They  calculate 
shear  stresses  using  equations  by  Hutchinson  et  al ■  and  Norton's  law,  and  they  introduce  Kachanov's  damage 
parameter  concept  by  writing  the  damage  accumulation  rate  as 


The  crack  is  assumed  to  grow  in  a  stepwise  manner  over  a  length  p  each  time  the  accumulated  damage  at  the 
crack  tip  reaches  a  critical  value. 

The  final  result  depends  strongly  on  the  relative  values  of  the  exponent  b  and  n  the  exponent  in  Norton's 
law. 


b/(n  +  1)  >  1 

-  n  b 

da  ,  f  2n  T0  )n+1 

dt  ’  A\it  (n  +  0  '  Tq”~ j 

Here  c  is  Riemann's  zeta  function. 
b/(n  +  1)  -  1 


<^t)- 


da  2n  T0  f(a  "  Vl  ,* 

dt  '  ir(n  +  1)  '  Y0  ’  n[  P  J  ' 


a  -  aQ  >>  p 


b/(n  ♦  1)  <  1 


'(  b  ) 

rf,  .  b  \ 

I  2n 

_^2_V 

U  ♦  t/  • 

r\  "  ♦  1/ 

lit(n  +  1 

T  y»  J 

(11.10) 
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Here  r  is  Che  gaoma  function. 

Freeman  (78)  calculates  Che  crack  propagation  rate  as  a  function  of  the  reference  stress  defined  as 

P 

0  m  ‘ 

ref  pB  w 
s 

where  u  is  the  usual  ratio  of  yield  loads  and  where  B  and  w  are  specimen  thickness  and  width. 

For  the  single  edge  crack  tension  specimen  (SENT) 


(11.11) 


('-=)! 


Experimentally  it  was  found  that  °ref  varies  with  time  according  to  the  equation 
°refP(tRO  '  C)  "  Z 


(11.12) 


(11.13) 


Here  tR^  is  the  rupture  time  associated  with  the  initial  value  of  o  ^  and  is  related  to  rupture  stress  oR 
according  to  the  empirical  formula: 

aRe.tR-Z  (11.14) 


The  result  is 


da  _  1  w0,5  .  P0-5  _  (2p-1)/2 

dt  "  2p  •  0.5Z  1  °ref 

s 

For  the  compact  tension  specimen  there  is  found  similarly: 

da  1  w0,5  .  P°'5  (2p-1)/2 

dt  1 . 15p  ’  0.5,  *  °ref 

D  L 


(11.15) 


(11.16) 


Using  concepts  from  the  reference  stress  method.  Freeman  relates  creep  and  elastic  CODs  of  the  SENT 
specimen  to  creep  and  elastic  strains  of  a  creep  specimen. 

Using  the  elastic  analogue  and  Norton's  law  he  then  calculates  the  creep  COD  rate  as  a  function  of 
reference  stress 


v  -  P  .  C(a/w) 
c 


^N  n-1 

B  '  °ref 
s 


(11.17) 


Here  C(a/w)  is  the  non-dimensional  compliance  and  BR  is 


the  Norton  law  constant 


C(a/w)  -  CEB  =  X  ,  EB 
S  P  s 


(11.18) 


Crack  propagation  rate  can  now  be  related  to  creep  COD  rate  for  the  SENT  specimen  as 


0.5„0.5  r  B  vc 


da  J_  w“'"P“‘"  [  s 
dt  “  2p  B  0.5Z  [bn  .  P  .  C(a 7 


lwTj 


2p-1 

2n-2 


(11.19) 


and  for  the  CT  specimen  as: 


da  _L_  w°-5P°-5 

dt  *  1-1 5p  •  g  0.5Z 


B  v  n 

s  c 

,Bn  •  P  •  C(a/w) J 


1 £li 

2n-2 


(11.20) 


Radhakrishnan  et  al ■  (46)  treat  the  single  edge  crack  bend  specimen  and  relate  the  incremental  COD  to 
the  increment  of  the  non-dimensional  crack  length  A  ■  (a/w)  as 


dd  -  Ak  .  a  .  m  .  d\  (11.21) 

tip 

where  A  and  k  are  constants. 

They  than  calculate  stresses  near  the  crack  tip  by  taking  four  terms  instead  of  the  usual  first  one  of  the 
LEFM  series  expansion  of  stress  as  a  function  of  distance  to  the  crack  tip  and  the  distance  from  the  crack 
tip  to  the  neutral  axis.  Bending  moment  and  load  P  are  introduced  through  the  equilibrium  equation.  This 
leads  to  the  equation 


d6 .  xg_(,  ♦  »)(4-i) 


2n  2n 

(GF)>  (dA) 


(11 .22) 


Here  n  is  the  exponent  from  Norton's  law  and  (GF)  is  an  expression  containing  gamma  functions  of  arguments 
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containing  n. 

The  authors  state  that  consequently  the  following  equation  must  also  hold 


and  that  finally 


2n 

m  2n-m 


(11.23) 


(11.24) 


Here  a  is  the  initial  value  of  the  nominal  bending  stress  due  to  P.  These  equations  are  the  basis  of  the 
correlSting  parameter  used  in  figure  17. 

Rice  et  al  (79)  have  analysed  stresses  and  strains  at  the  tip  of  a  growing  crack  on  the  basis  of 
Prandtl's  slip  line  field  theory.  They  present  the  formula 


d  .  a  d  0  da  .  R 

dt  r  ”  Cg  dt  E  dt  n  r 


(11.25) 


Here  6  is  the  COD  at  a  distance  r  behind  the  crack  tip.  The  constants  a,  6  and  R  were  determined  by 
comparison  to  results  of  a  FEM  analysis.  It  was  found  that 

a  *  0,65,  6  *  5,08  and 


R  -  0,2 


(11.26) 


Purushothaman  and  Tien  (80)  assume  an  Orowan-Irvin  elastic  stress  distribution  at  the  crack  tip.  They 
take  it  that  for  any  line  element  dx  in  front  of  the  crack  tip  there  will  be  a  rupture  time  tR(x) 
determined  by  the  local  stress  and  that  the  crack  propagation  rate  can  be  calculated  as 


da  _  dx 

dt  dtR  ^at  the  crack  tip 


(11.27) 


Rupture  time  is  calculated  from  the  empirical  equation 


*R  '  VW 


(11.28) 


and  minimum  creep  rate  from  the  phenomenological  equation 


•  n  fo(*)\n  "c 

emin  “  BN|E(T)J  eXp  RT 


(11.29) 


The  result  is 


da  ^SL  I"  “f  K  V™  -°Qc 
dt  "  anB_  '  P  lvp  0.5/  exp  RT 


(11.30) 


Here  p  is  the  crack  tip  radius  and  Y  is  the  function  Y(a/w)  from  the  relation 


K  ■  Y  .  o(ira) 


(11.31) 


It  should  be  noted  that  the  use  of  the  debatable  relation  between  rupture  time  and  minimum  creep  rate  is 
not  necessary  since  one  can  use  a  phenomenological  equation  directly  relating  rupture  time  to  stress  and 
temperature. 


Van  Leeuwen  (2,  3)  has  considered  the  diffusion  of  vacancies  towards  the  crack  tip  and  related  the 
chemical  potential  of  a  vacancy  at  the  crack  tip  to  the  rate  of  energy  dissipation  G. 

The  result  is 


—  ■  —  (C  -  C  ) 
dt  R  (Ct  w' 


(11.32) 


Here  (Ct  -  Cv)  is  the  difference  between  the  vacancy  concentrations  at  the  crack  tip  and  the  crack  wall 


C  -  C 
t  w 


-<*i) 


+  G(a)0 


(11.33) 


Cq  is  the  vacancy  concentration  at  infinity 

puis  numerically  equal  to  the  nominal  tensile  stress  normal  to  a  circular  embedded  crack  with  diameter  2a. 
Y  is  the  surface  tension 
p  is  the  crack  tip  radius 
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fl  is  Che  atomic  volume 
D  is  the  diffusion  constant 

R  is  the  average  distance  between  crack  nuclei 
k  is  Boltzmann' 8  constant 
T  is  the  absolute  temperature 

Jq  and  J  are  Bessel  functions  of  the  first  kind 
is  the  first  root  of  the  equation  J((X)  ■  0 

Henshall  et  al.  (31)  consider  creep  crack  growth  by  vacancy  diffusion  using  equations  derived  by  Vitek 
and  by  Speight  and  Beere  respectively. 

In  the  first  case  there  is  found 


(11.34) 


(11.35) 

Dg  is  the  coefficient  of  self  diffusion  for  the  grain  boundary 
6°  is  the  width  of  the  grain  boundary 
Si  is  the  atomic  volume 
u  is  the  shear  modulus 
k  is  Boltzmann's  constant 
T  is  the  absolute  temperature 

d  is  the  width  of  the  crack 

c 

In  the  second  case  there  is  found 


da  _  0.S16  X  (KY 
dt  *  d  4  W 
c 

Here 


(11.36) 


Here  n  and  B  are  constants  from  Norton's  law  and  y  is  the  surface  tension. 


12.  DISCUSSION,  CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  found  that  linear  elastic  fracture  mechanics  can  be  applied  sucessfully  to  the  residual 
strength  of  brittle  materials,  crack  growth  due  to  stress  corrosion  and  crack  growth  due  to  fatigue.  Even 
then  the  stress  intensity  factor  K  is  not  the  sole  parameter  determining  residual  strength  or  crack 
propagation  rate.  The  relationships  are  influenced  by  some  additional  geometric  and  environmental  para¬ 
meters.  Therefore  it  is  not  realistic  to  expect  that  the  growth  rate  of  creep  cracks  will  be  determined  by 
one  single  fracture  mechanics  parameter. 

With  the  methods  to  account  for  plasticity  in  residual  strength  problems  a  return  is  discernible  from  the 
stress-oriented  K-based  approached  to  the  energy-oriented  approach  based  on  G,  R  or  J.  Therefore  it  may  be 
expected  that  the  crack  growth  rate  under  creep  conditions  will  correlate  well  with  a  parameter  directly 
related  to  the  dissipation  of  energy. 


In  the  case  that  a  stress-oriented  approach  is  selected  nevertheless,  the  redistribution  of  stresses 
due  to  creep  and  stress  relaxation  has  to  be  accounted  for.  The  K  factor  then  loses  the  greater  part  of 
its  relevance. 

Notable  is  the  erratic  success  of  K,  onet  and  vrej  as  a  correlating  parameter  for  crack  growth  rate. 

The  expectation  is  not  borne  out  that  in  the  cases  of  nominally  brittle  materials  or  narrow  intercrystalline 
cracks  the  K-factor  will  do  better  than  cnet  or  °ref'  sPec^men  geometry  plays  an  important  role.  If  it 

is  such  that  lateral  contraction  at  the  crack  tip  is  highly  constrained,  then  the  energy  dissipated  will 
be  spent  much  more  on  crack  propagation  than  on  plastic  deformation  and  the  crack  propagation  rate  will 
correlate  fairly  well  with  K.  If  correlation  is  bad,  this  means  that  the  relation  of  da/dt  with  K  or  °net 
is  influenced  heavily  by  specimen  shape  or  dimensions  and  the  initial  values  of  crack  length  or  load. 
Sometimes  a  better  correlation  is  obtained  by  working  with  the  relative  crack  length  a/w  or  by  incorpo¬ 
rating  rupture  time  or  load  in  the  correlation  parameter.  These  seem  to  be  tricks  that  will  work  only  under 
limited  conditions. 

Because  COD  and  COD  rate  are  closely  related  to  energy  dissipation  it  can  be  expected  that  they  will 
correlate  we 1 1  with  the  crack  propagation  rate,  it  being  understood  the  re-initiation  of  a  crack  will 
probably  occur  at  a  critical  value  of  5  and  that  afterwards  da/dt  will  correlate  with  dS/dt.  In  practice 
this  will  be  jf  little  use  because  of  the  lack  of  fast  methods  to  calculate  6  and  d6/dt  for  a  crack  in  a 


structun 1  element,  with  the  aim  to  predict  da/dt.  Some  use  can  be  made  nevertheless  of  the  approximate 
expressions  for  G,  J  and  J*  applicable  to  certain  specimens  and  containing  6  and  d£/dt. 

A  static  quantity  like  K,  o  .,  G  or  J  seems  to  be  better  suited  for  predictive  purposes  than  a  dynamic 
quantity  such  as  d6/dt  or  J*. 

Coleman  et  al.  (81)  have  been  fairly  successful  in  correlating  the  growth  rate  of  creep  cracks  in  pressure 

vessels  with  K  and  a  . 

It  is  remarkable  thaErelatively  few  attempts  have  been  made  to  correlate  da/dt  with  J,  and  that  much  more 
is  expected  from  the  dynamic  parameter  J*.  It  should  be  considered  that  in  order  to  calculat'  J*  one  needs 
values  of  creep  rates  and  deformation  rates  and  that  these  in  their  turn  will  be  determined  by  static 
quantities  like  a  ,  or  J.  Then,  nobody  challenges  the  fact  that  creep  rate  will  be  determined  by  a  static 

quantity  like  strisiS  a. 

It  is  reconmended  therefore  to  put  more  effort  in  attempts  to  correlate  da/dt  with  J. 


An  interesting  development  is  that  attempts  are  already  made  to  correlate  the  crack  growth  rate  in  high 
temperature  fatigue  with  J.  (..  \dananda  and  Shahinian,  82).  It  is  a  fortunate  fact  that  so  many  alternative 
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and  fast  methods  are  available  to  calculate  J  and  J*.  The  extent  of  agreement  in  their  results  is 
encouraging. 

In  addition  it  is  felt  that  parameters  like  dj/da,  d6/da  and  even  dj*/da  merit  more  attention. 

The  literature  has  shown  that  da/dt  correlates  rather  well  with  J*  although  the  specemen  geometry  still 
seems  to  have  some  influence. 

Riedel  (83)  and  Riedel  and  Rice  (84)  define  a  creep  zone  in  analogy  with  the  plastic  zone.  Its  boundary  is 
the  locus  of  points  where  creep  strain  equals  elastic  strain.  The  creep  zone  grows  with  time. 

Riedel  and  Rice  state  that  at  the  start  of  the  creep  cracking  process  the  creep  zone  will  be  small  and  that 
then  K  will  be  determining.  Lateron  the  creep  zone  will  be  relatively  large  and  then  J*  will  be  determining. 
The  present  author  is  rather  opposed  to  a  discontinuous  approach  like  that.  It  is  the  demonstration  of  a 
desire  to  keep  on  using  the  K-factor  whenever  possible  and  pushing  it  to  or  even  beyond  its  limits  of 
applicability.  It  should  be  remembered  that  the  K-factor  strictly  speaking  is  relevant  only  under  the 
conditions  that 

-  the  deformations  are  elastic 

-  the  crack  is  not  too  small  relative  to  typical  dimensions  of  the  cracked  body 

-  the  problem  is  two-dimensional. 

If  these  conditions  are  not  fulfilled  it  is  better  to  depart  from  the  stress-oriented  K-based  approach  and 
to  return  to  an  energy-oriented  approach  based  on  parameters  like  J  or  J*. 

If  the  conditions  are  fulfilled  one  can  nevertheless  take  the  energy-oriented  aporoach  since  then  K  is 
so  closely  related  to  G  and  J. 

A  very  important  complication  in  the  application  of  fracture  mechanics  to  creep  cracks  is  that  the 
materials  properties  of  the  structures  considered  are  influenced  so  heavily  by  processing  (welding  and  heat 
treatment)  and  by  accumulated  creep  and  high-temperature  exposure  (cavitation  and  ageing).  Therefore,  if 
creep  growth  parameters  are  determined  experimentally,  this  will  have  to  be  done  on  material  that  has 
undergone  a  representative  pre-treatment.  Examples  can  be  found  in  Sivems  and  Price  (42,  85), 

Neate  and  Sivems  (20),  Neate  (22),  Gooch  (26,  27,  48),  Gooch  and  King  (86)  ,  Ritter  and  Formby  (87), 

Worswlck  and  PiTkington  (41,  52,  88),  Christian  et  al ■  (59) .  ~ 

Finally  it  has  to  be  noted  that  in  the  present  review  of  the  literature  it  has  been  assumed  tacitly 
that  failure  of  a  high-temperature  structure  is  initiated  by  formation  and  growth  of  a  single  or  dominant 
crack.  This  will  not  always  be  the  truth.  Another  failure  mode  consists  of  creep  cavitation,  i.e.  the 
formation  of  micro-pores  at  critical  locations.  These  voids  will  grow  by  plastic  deformation  and/or 
vacancy  diffusion  until  the  cross  section  has  been  weakened  to  such  an  extent  that  rupture  occurs  by  over¬ 
load.  A  critical  review  of  the  literature  concerning  void  growth  by  vacancy  diffusion  has  been  given  by 

Van  Leeuwen  among  others  (2) .  A  model  to  predict  the  residual  strength  of  a  material  weakened  by  voids 
has  been  developed  by  Scott  et  al.  (89). 

Relevant  to  this  problem  is  the  analysis  by  Van  Elst  (90)  who  considered  an  array  of  defects  or  cracks 
rather  than  a  single  or  dominant  one.  He  has  based  his  approach  on  an  analysis  by  Koiter  concerning  a 
linear  periodic  array  of  through  cracks  of  size  a  and  distance  d  leading  to  a  representatitve  stress 
intensity  factor 

K  -  o(2d  .  tan  ~  A)0,5  (12.1) 

where  A  •  (a/d).  A  finite  element  analysis  of  a  two  dimensional  array  of  parallel  through  crack,  length  a, 
at  mutual  distances  d.  »  d^  a  d  has  shown  a  similar  tangent-like  dependence  of  the  effective  K  factor  on 
A  ”  (a/d).  This  analysis  is  now  being  extended  to  three-dimensional  arrays. 
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TABLE  2 

Results  of  the  calculations  by  Goldman  and  Hutchinson 
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Fig.  24  Test  results  by  Ohji  for  type  304  stainless  steel  at  650  C  (4). 
N.B.  dashed  curves  pertain  to  CT  specimens 
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Fig.  25  Test  results  by  Ohji,  Ogura,  Kubo  and  Katada  (57) 
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Fig.  26  Correlation  of  da/dt  with  0net  and  J*  for  tests  on  CNC-M  specimens  at  various  stress  levels 
(upper  plots),  on  CNC  and  RNB  specimens  of  different  sizes  at  various  stress  levels  (middle 
plots),  on  CNP  specimens  of  several  sizes  at  one  stress  level  (lower  plots).  Taira  et  al.  (58) 
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showing  effects  of  primary  creep 
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